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We investigate a model system for the injection of fermionic particles from filled source sites into an empty 
chain. We study the ensuing dynamics for Hermitian as well as for non-Hermitian time evolution where 
the particles cannot return to the bath sites (quantum ratchet). A non-homogeneous hybridization between 
bath and chain sites permits transient currents in the chain. Non-interacting particles show decoherence in 
the thermodynamic limit: the average particle number and the average current density in the chain become 
stationary for long times, whereas the single-particle density matrix displays large fluctuations around its 
mean value. Using the numerical time-dependent density-matrix renormalization group (t-DMRG) method 
we demonstrate, on the other hand, that sizable density-density interactions between the particles introduce 
relaxation which is by orders of magnitudes faster than the decoherence processes. 
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1 Introduction 

Newton's equations describe the deterministic time evolution of a given initial state of classical particles. 
In principle, the particle velocities and positions are known for all times. In practice, in isolated systems 
and for long times, macroscopic observables which involve the averages over many particles can be calcu- 
lated equally accurately starting from thermodynamic considerations. The 'thermalization' of macroscopic 
systems is one of the building principles of statistical mechanics. According to statistical mechanics, ex- 
pectation values for macroscopic observables should, for long times, become independent of time and, 
moreover, become equal to their statistical average; for a recent example of the relaxation of classical par- 
ticles in a one-dimensional box, see Ref. |1|. The 'thermalization principle' equally applies to quantum 
mechanics where the Schrodinger equation describes the deterministic evolution of a given initial state in 
time. Recently, the question under which conditions thermalization is possible in closed quantum systems 
has also been studied, both experimentally and theoretically |2-4|. According to the eigenstate thermal- 
ization hypothesis, each generic state of a closed quantum system already contains a thermal state which 
is revealed during unitary time evolution due to dephasing and relaxation processes. However, proofs for 
this hypothesis have so far only been obtained in very specific settings l5]|6|; its appUcability for generic 
quantum systems remains a matter of current research. 

In this work, we are interested in the quantum-mechanical injection of particles into a chain. For this 
toy model we can study the approach to a steady state analytically and numerically. The time evolution 
of prepared initial states in one-dimensional quantum systems such as our toy model can be studied ex- 
perimentally using cold atomic gases on optical lattices confined to one dimensional potentials |l2][3]|7][8l . 
These systems are to a very high degree isolated from the environment and are therefore ideally suited to 
study non-equilibrium dynamics in closed quantum systems. Besides the Hermitian time evolution we also 
address the non-Hermitian case for our toy model where the chain particles cannot return to the bath sites 
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(a 'quantum ratchet' |l9][T0l). In general, effective non-Hermitian Hamiltonians are naturally obtained if 
one simulates a dissipative environment for a quantum model in the form of a Lindblad equation ifTTI . An 
effective non-Hermitian Hamiltonian has also been used to study the depinning of flux lines in supercon- 
ductors lfT2l . Interesting from a fundamental point of view are, in particular, so-called pseudo-Hermitian 
Hamiltonians — including Hamiltonians having PT-symmetry — which have real spectra and therefore of- 
fer a natural framework to extend standard quantum mechanics IfTTI jTS I. 

In Sect. |2] we recapitulate the Schrodinger equation for time evolution in quantum mechanics. We treat 
both cases of Hermitian and non-Hermitian Hamilton operators. In Sect.[3]we introduce our model system 
for spinless fermions. We consider a chain of sites which is initially empty and is coupled to initially filled 
bath sites. At finite times, the hybridization introduces particles in the chain where they can move between 
neighboring sites. The hybridizations can be homogenous or inhomogeneous; in the latter case they induce 
a transient current in the chain. 

In Sect. |4] we discuss the particle density, the single-particle density matrix, and the current density for 
non-interacting fermions. In the thermodynamic limit, the particle and current densities become constant 
for large times because the contributions of individual modes interfere destructively ('decoherence'). The 
oscillations in the single-particle density matrix on the other hand reveal the absence of relaxation for 
non-interacting particles. Furthermore, we investigate the influence of the ratchet condition on the chain 
dynamics. In Sect. |5]we include a density-density interaction between the particles. We use the time- 
dependent density-matrix renormalization group (i-DMRG) method to evolve the quantum-mechanical 
states in time. The interaction between the particles induces a fast relaxation to stationary distributions for 
the particle density, single-particle density matrix, and current. For our model system, with interactions of 
the order of the bandwidth, we find that the relaxation time scale is two orders of magnitude shorter than 
the time scale for dephasing. 

A summary and conclusions. Sect. |6] , close our presentation. Calculations for the asymmetric chain 
filling in the ratchet case are deferred to appendix lAl 

2 Hermitian and non-Hermitian time evolution 

We begin our discussion with the Schrodinger equation which defines the time evolution for Hermitian 
and non-Hermitian Hamilton operators. As in standard quantum mechanics, expectation values are to be 
taken using the time-dependent wave function. The von-Neumann equation of motion of the corresponding 
statistical operator can be cast into the general form of a Kossakowski-Lindblad equation. 

2. 1 Schrodinger equation 

We consider a finite-dimensional Hilbert space V. Its elements are the states \ip). We use the standard 
scalar product {4>\tp) S C where {(j>\ is the shorthand notation for the linear functionals {(f)\-) which form 
the dual space V* of our Hilbert space. For an operator O which acts in V, its adjoint is defined by 
{4)\d^ IV') = (V'|0|V))* for afl g V* and IV^) e V where z* is the complex conjugate of z S C. 
For time t > 0, the time evolution of a state |?/') and of its adjoint state (?/'| is given by 

mt)) = exp (-iiJt) 1^) , {m\ = ('^l exp (iiJ^t) , (2.1) 

where _ff is a time-independent Hamilton operator and W is its adjoint. 

In the following we shall assume that the spectrum of H is real, i.e., the eigenstates of H obey 

H\En)=E„\En) (n=l,2,...) (2.2) 

with En G M. As shown in Ref. |13| this implies that H is pseudo-Hermitian and H^^ has the same 
spectrum, 

H^\E;;)^En\K) (n = l,2,...). (2.3) 
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Eq. ( 12.31 1 implies that there is a left eigenvector {En \ of H for every right eigenvector of H with the 
same eigenvalue £"„. 

Note that {E^\ = only if H is Hermitian, = H. The eigenstates of H and H^i for different 

energies are orthogonal {En\E„i) = if £^„ ^ Em- Using proper linear combinations of eigenstates for 
degenerate eigenvalues we can write for the Hamiltonian, its adjoint, and the unit operator 

= ^£;„|s„)(:b:| , i/^ = ^^„|:e;)(£;„| , i = ^ = ^ |^„)(:b:| = i^ 

71 n n n 

(2.4) 

where (i?„|i?„) = 1 thus forming a biorthonormal set of eigenvectors. The Hamiltonian and its adjoint 
obey the (pseudo-)Hermiticity condition lfT3l 

fjH^H^f, , fi = Y,\E:){E^\^V^ , 77-1 -^|S„)(^«| . (2.5) 

n n 

Note that all definitions reduce to the standard expressions in the case of Hermitian time evolution with 

H = wAfj^ i. 

2.2 Measurements 

In standard quantum mechanics, the expectation value for an operator O is given by 

,o),t)^m9\m (2 6) 

^^^"^'^ {m\m) • ^ ^ ^ 

For pseudo-Hermitian time evolution, this definition remains unchanged. It is important to note that the 
time evolution described by ( 12. Il l is not unitary so that the norm {il}{t)\il}{t)) is not conserved and must be 
considered separately. 

As a consequence, the equation of motion for the statistical operator changes its form. The statistical 
operator for the pure state \'ip{t)) is defined by 

m))m)\ ,^ 7, 

With the help of the statistical operator, expectation values can be expressed in the form 

(6)(t) -Tr(p(i)d) . (2.8) 

where Tr denotes the trace which implies the sum over the expectation values of an orthonormal basis set. 
The statistical operator obeys — p, (? — p and Tip = 1. Its time dependence is given by the equation 
of motion (generahzed von Neumann equation) 

i^p(i) = ilp{t) - p{t)H^ + p{t){H^ - H){t) . (2.9) 

When we split the Hamiltonian into a Hermitian part = Hq and a non-Hermitian perturbation V ^ V'^ 
we can equivalently write 

^p(t) ^ -i [Ho,m] _-i{{v- {v)it))m miv^ - (v^m)) • (2.10) 

Apparently, the generalized von-Neumann equation takes the form of the Kossakowski-Lindblad equa- 
tion l,16.,17J . 

Copyright line will be provided by the publisher 



4 



F. Gebhard, K. zu Munster, J. Ren, N. Sedlmayr, J. Sirker, and B. Ziebarth: Particle injection into a chain 



3 Particle injection into a chain 

As an example we study the injection of fermionic particles from filled source sites into an empty chain. 
We start with a translationally invariant coupling of the source and bath sites. Next, we consider an in- 
homogeneous coupling which permits us to generate a transient current density in the chain. For the case 
of non-interacting fermions, analytical expressions can be worked out for the particle density, the single- 
particle density matrix, and the current density. For interacting fermions, we evolve the states in time with 
the density-matrix renormaUzation group (t-DMRG) method. 

3 . 1 Hamilton operator 

We consider a finite chain with L sites on which spinless fermions can move between neighboring sites, 

L-l 

T = -J^ (4+iQ+c|q+i) , (3.1) 

where c\ (ci) creates (annihilates) a fermion at site I. Later, we will set J = 1 as our energy unit (bare 
bandwidth W = 4). 

The chain fermions are locally hybridized with reservoir fermions, 

L 

1=1 

The hybridization strength for the fermion transfer from the reservoir to the chain, 7/ in, can be different 
from the transfer from the chain to the reservoir, 7z,out. In this work we focus on the Hermitian case, 
7;. in = 7/, out, and the ratchet case, 7/^out = 0, where the chain fermions cannot return to the reservoir sites. 
We will study the translationally invariant case, 7/ in = 7, and the inhomogeneous case with 7/+i,in < 
liM = J exp(-;/L). 

In Sect. |5] we consider numerically the time evolution in the presence of density-density interactions 
between the fermions, 

y = ^ 14 (c|q - 1) {cl,ci+, -l)+Vs i4s^ - i) {sl,§i+^ - i) + T/e. (4q - i) {4 k -I). (3.3) 
I 

All interaction parameters are chosen to be positive. The total Hamiltonian is given by 

H = f + S + V . (3.4) 

The model is shown pictorially in Fig.[T] 

f f 

171 1 72 

6 6 

J 

Fig. 1 Model for the particle injection into a chain. The solid circles represent the source sites, the open circles the 
chain sites. The solid and dashed lines denote the bonds with hopping amplitudes J and 7; respectively. 
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3.2 Initial state and measured quantities 
In order to study the particle injection from the reservoir into the chain, we choose our initial state as 



|V(t = 0)) = l^o) =n^;Vac) . (3.5) 



1=1 



At time t = all source sites are occupied and the chain is empty. 

We are interested in the average particle density in the chain, ndt), the single-particle density matrix of 
the chain fermions, Ap^^it), and the current density, jc{t)- These quantities are defined by 



ApAt) = {alak)it) , flfe- y^T^I^sinf-^— (fc = 1, 2, . . . , L) , 



1=1 

L 



ricit) = (3.6) 

k=l 

i ^"^ 2i ^ 

jc{t) = 2 ^{cl+iCi - c|c/+i)(t) = -J—^ i(^'P)^p,fcW ' 

1 = 1 p.k=l 

p^k 

(1 - (-l)P-^) sin(p7r/(£ + 1)) sin(fc7r/(L + 1)) 



j{k,p) = 



i{kTT/{L + 1)) - cos(p7r/(L + 1)) 



For a homogeneous coupling, ji — 7, the single-particle density matrix is diagonal, ^(i) = Sk,pAk,k, 
and the current density is zero for all times, jc{t) = 0. 

The long-time mean of a time-dependent quantity X{t) is defined as 

X{t) = lim - [ X{t)dt . (3.7) 

T^OO T Jq 



4 Non-interacting fermions 

First, we address the question of how effective the injection of non-interacting fermions into the empty 
chain is. To this end, we consider the average particle number, the single-particle density matrix, and the 
current density in the chain, as functions of time. 

At short times we expect a linear increase of the particle density in the chain because we start with 
all particles at the source sites. At later times, we shall see an oscillatory behavior of the particle density 
because of the coherent quantum-mechanical movement of the particles within and into the chain (and back 
onto the source sites for the Hermitian case). In the non-interacting case, a constant average occupation is 
reached only for long times because the two-level systems stay coherent for relatively long times so that 
the oscillations in the particle density die out only slowly. Interactions dampen these oscillations strongly 
and also modify the value in the long-time limit considerably, see Sect. |5] 



4. 1 Homogeneous chain filling 

We begin with the translationally invariant case for non-interacting particles, ji^in — 7in and 7;_out — Tout 
in eq. (13.2b . For a homogeneous chain filling, the single-particle density matrix is diagonal and the current 
density is zero for all times. 
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4.1.1 Hermitian time evolution 
In the Hermitian time evolution with 7i„ = 7out = 7 the model reads 

L-l L 

i=i j=i fc 

where bk denotes the Fourier transform of the operators s; defined analogously to the Fourier transform of 
the chain operators in eq. ( 13.61 1. The dispersion is given by et = —2 J cos k with fc = 7rn/(L + 1). 

Heisenberg's equations of motion provide the simplest way to solve this problem analytically. The 
resulting coupled differential equations read 

CLk = i[H, ak] = -i{skak + jh) and 

bk = i[H, bk] = -h&k . (4.2) 
It is straightforward to write down a formal solution. 



where 



bk{t) 



^'^^^^^ =e--^*/2Mk(0('?) (4.5) 



Re-expressing the exponential matrix leads to 

'ak{t)\ ^ 

Mt)J ^ Kbk 

with 

Mk.i2^Mk,2i =- -i^sm{Ekt/2) , (4.6) 

Ek = 

The zero-time correlation functions are very simple, 

(V^ol4afcl^o) = (V-olalA-l^o) =0 , (Vol&I^felV^o) = 1 ■ (4.7) 
The single-particle density matrix is diagonal and given by 

Ak.kit) = \Mk,i2f = % sin^ {Ekt/2) . (4.8) 

^k 

It describes the beating of particles in the two-level systems with energies ek-i,2 — (l/2)(£fc ± Ek) which 
result from the hybridization of the bath particles at zero energy with the chain particles with dispersion 
£k- 

The average particle number in the chain becomes 

n,{t) = jJ2%^^^'iEkt/2). (4.9) 



El 
k k 
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Its long-time mean is straightforwardly extracted, 



2L^ El 



(4.10) 



In principle, the oscillations around the mean will not die out for long times. For any finite chain there exists 
a recurrence time t^^^ ^ such that for the numerator of the sum ( 14.9b the relation sm^{tl.^^Ek/2) < t/L 
holds for all allowed /c-values. However, this Poincare time is very large even for systems L — because 
the energies Ek are not commensurate. 

In the thermodynamic limit (L oo) eq. ( I4.9l l reduces to 



2 

ndt) = - / 

Jo 



de 



1 



(4.11) 



In this case the oscillations in nc{t) do vanish for t —?' oo because the infinitely many two-level systems 
dephase completely. We find 



ricit — >■ oo) 



1 



1 



2 VI + (Jhy 



(4.12) 



For J/7 — > only half of the chain sites are occupied on average. Recall that, for J = 0, we describe a 
collection of independent two-level systems where each particle oscillates between two levels such that it 
is equally likely to be found on the chain and on the bath site. 
For J/7 <C 1, the integral (I4.11l i can be solved to leading order. 



n,it) « i - i Jo(t JV(27)) [cos(27i) cos(aV(27)) + sin(27t) sin(aV(27))] . (4.13) 

Here, the trigonometric functions describe the fast oscillations around the mean value whereas the Bessel 
function Jo(tJ^/27) gives the envelope. The oscillations around the mean value decay in the long time 
hmit proportional to y/2j/JH. The slow oscillatory decay proportional to 1/ V< results from the square- 
root divergence in the single-particle density of states in one dimension. 



0.08 
0.06 
\ 0.04 
0.02 

°, 
0.6 

n 0.4 

c 

0.2 


r' 
0.! 

n 0.( 

0.1 

0, 



' 1 


1 


' 1 












■/ , 1 


1 


1 


1 


10 


20 


30 


40 5( 




10 


20 


30 


40 SO 


T 


1,1.- 



0.8 



0.6 



0.4 



0.2 



(d) 



1 



/ V/ - 







50 



100 200 300 400 500 

Jt 



Fig. 2 Particle density in the chain nc(t) for non-interacting fermions (Hermitian time evolution). Left: a chain with 
L = 50 sites and (a) 7/ J = 0.1, (b) 7/J = 1, (c) 7/J = 10; (d): thermodynamic limit for 7/J = 10 (circles: 
time steps &t — 0.1). The solid line is the long-time limit, the dashed lines show the envelope function according to 
eq. l l4l3> . 

In Fig. 12] we show results for different ratios of 7/ J for a chain with L = 50 sites for times Jt < 50. 
For a finite chain, the average particle number will never become constant. Instead, revival oscillations 
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will occur because the Hilbert space has a finite dimension. In the thermodynamic limit, the oscillations 
slowly decay to a constant average particle number in the chain given by eq. ( I4.11l i. see Fig. EJd), plotted 
for 7/ J = 10 and times Jt < 500. 

As a comparison, we show results for the ratchet case in Fig.[3]for the same values of = 7 (7out = 0) 
and different time spans. Even in the ratchet case, not all particles are injected into the chain at long times. 
Only for 71,1 — > 00 does the particle number in the chain approach unity, see eq. ( 14.191 ). In addition, one still 
observes fluctuations so that the particle density in the chain is not a monotonically increasing function. The 
motion of the fermions on the chain leads to Pauli blocking and to quantum-mechanical interference effects 
which result in the observed incomplete and oscillatory chain filling. In comparison with the Hermitian 
time evolution it can be seen that in the ratchet case the quantum-mechanical fluctuations are smaller and 
the chain filling is bigger for the same parameter 7/ J. 
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Fig. 3 Particle density in the chain nc{t) for non-interacting fermions on a chain with L = 50 sites for (a) 7in/ J = 
0.1, (b) 7in/ J — 1, and (c) "fin/ J = 10 (ratchet time evolution). The straight lines denote the mean values. Note the 
different scales on the axis. 



In Fig.|4]we show the single-particle density matrix Ak^k for Hermitian time evolution with -f/J—l 
at various times. At small times all particles enter the chain at roughly the same moment which leads 
to a relatively flat distribution for Jt = 1. At larger times, the individual dynamics of each fc-point 
becomes visible and a lot of maxima and minima have developed already by Jt ~ 50. At very large times, 
Jt = 5000, the map Ak^k{t), corresponding to a large collection of undamped oscillators, displays a non- 
monotonic ("chaotic") distribution. The size of the fluctuations around the time-averaged distribution is as 
large as its mean. The system does not relax. 

The corresponding results for the single-particle density matrix Ak^k for non-interacting electrons in the 
ratchet case are shown in Fig.|5] The average distribution is sharply peaked at fc = 7r/2, see eq. (14.18b . 
This means that for large times the fc-points near 7r/2 show very little variation with time while the fc-points 
away from 7r/2 are still strongly oscillating in time. 
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Fig. 4 Density matrix Ak,k{t) in the Hermitian time evolution with 7/ J = 1 for a chain of length L — 1000 as a 
function of the discrete momenta k = 7rn/(L + 1) for times Jt = 1, 10, 50 (left) and Jt — 5000 (right) where the 
solid line is the long-time mean. 




Fig. 5 Density matrix Ak.k (t) in the ratchet case with 7in/ J = 1 for a chain of length L — 1000 as a function of the 
discrete momenta k = nn/ {L + 1) for times Jt — 1, 10, 50 (left) and Jt = 5000 (right) where the solid line is the 
long-time mean. 



4. 1 .2 Ratchet case 

Here we provide the analytical expressions on which the data in Figs. [3] and |5] are based. We treat the 
translationally invariant ratchet case, 7/, out = and 7;_in — 7. The Fourier transformed Hamiltonian reads 

H = J2^Skalak + -/alh] = Ho + S , (4.14) 

k 

where Hq is the Hermitian and S is the non-Hermitian part. In the interaction picture representation we 
find 

^(i) -7^e'^'=*4Sfe . (4.15) 
fc 

In this case one can explicitly calculate the ^-matrix because S{t) commutes with itself at different times 
and therefore one can directly integrate over times in S{t). The wave function becomes 

m)) = l[[^ + 19k{t)bkai]mO)) : 9k{t) = . (4.16) 

k=l 
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For the single-particle density matrix we find 

l + T^k-WP 4/(472) + sin2M/2) • ^ - ' 

In contrast to the Hermitian time evolution, the bare energy levels of the bath and the chain fermions do 
not hybridize. Therefore, the bare dispersion Ek appears as the beating frequency. Also note the time- 
dependence of the denominator which is absent for the Hermitian time evolution, see eq. (14.81 ). 

Again, we note that it is always possible to find a recurrence time tf.^^ such that the average particle 
number in the chain nc{t^c^) — J^k ^fc-fc(^roc)/-^ < ^ for ™y given e. At first, these recurrences seem 
counter-intuitive because, in the ratchet case, the amplitude for particles going back from the chain onto 
the source sites is zero. This clearly demonstrates the quantum character of the problem. 

The time averaged distribution function is given by 



= (4.18) 

which is sharply peaked at k ^ tt/2. In the thermodynamic limit we therefore find that the density of 
particles in the chain is given by 

nJt — ^ 00) = 1 — — arctan | — ) . (4.19) 

For a collection of independent two-level systems of bath and chain sites (J — 0), the chain is completely 
filled for large times because the particles cannot return to the bath sites. For J/7 7^ 0, however, we have 
nc{t — > 00) < 1, i.e., the chain is not completely filled despite the ratchet condition. Due to the motion of 
the particles on the chain, some bath particles are prevented from falling down into the chain because their 
chain sites are blocked by other particles (Pauli exclusion principle). 



4.2 Inhomogeneous chain filling 

For an inhomogeneous coupling, the problem must be treated numerically. It involves the diagonaUzation 
of the Hamiltonian (Hermitian time evolution) or the inversion of the self-energy matrix (ratchet case). 

4.2.1 Hermitian time evolution 
We rewrite the Hamilton operator as 

2L 

H=J2 H^ AP, (4-20) 

with P2m-1 = S„i, p2m = Cm for m = 1,2, ... ,i and H2m-l,2m — H2m,2m-1 ~ Im, H2m+2,2m — 

H2m,2m+2 = — J for TO = 1, 2, .... L — 1. The Hamiltonian is diagonal in a new single-particle basis, 

2L 2L 2L 2L 

Ar=l j=l N=l N=l 

where we used in the last step that the transformation matrix U, which is determined numerically, is real. 
The single-particle density matrix in momentum space is best calculated from the corresponding expression 
in position space. 




(4.22) 
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with 



ciAt) = {m\pliP2r\m) 



2L 



(En — EM)i 



2 m — 1, m — 1 , A-/ 



U2l,NU2r,M 



(4.23) 



Here we used the time evolution ( 12. Il l and the fact that all bath sites are occupied at time i = so that 
(?/'(0)|p2m_iP2n-i|V'(0)) = Sm,n- The Calculation of Ap,q{t) involves five lattice summations which are 
performed numerically. 

For an inhomogeneous coupling, ji = J exp{--l/ L), the particle-density in the chain ndt) is qualita- 
tively similar to the case of a symmetric filling for 7 = J. In Fig.|6] we show the filling of the chain for 
the unitary time evolution and the ratchet case for times Jt < 100. 
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Fig. 6 Particle density in the chain nc{t) for non-interacting fermions with hopping 7; — J exp{~l / L) on a chain 
with L = 200 (L — 1000) sites for unitary time evolution (left) and ratchet time evolution (right). 



In the unitary case, the fluctuations in the particle density are much smaller for the inhomogeneous 
hybridizations than for constant hybridizations. In addition, larger system sizes reduce the fluctuations 
further as can be seen from the comparison for system sizes L = 200 and L — 1000. In any case the 
particle density reaches its long-time value for times Jt > 1000 with small fluctuations around it. These 
reflect the transient current density in the system, see Fig.|8] In the ratchet case, we observe that the particle 
density in the long-time limit for the same set of parameters is higher than in the unitary case. Furthermore, 
the oscillations around this long-time mean are more regular and decay algebraically. 

In the case of an inhomogeneous coupling, 7; = Jexp(— Z/i), the single-particle density matrix is no 
longer diagonal. The solution of the equations (14.231 1 and ( |4.29t for L = 200 sites shows that dephasing 
is very efficient for the non-diagonal elements Ak,p^k- Even at short times, Jt ^ 1, the single-particle 
density matrix is almost diagonal with Ak^p=£k (t) being an order of magnitude smaller than Ak.k (t) for all 
k,p. Further decay in time of the off-diagonal elements is only algebraic thus leading to an algebraically 
decaying current density, see the following discussion. In Fig.|2]we show the time evolution of the diagonal 
components of the single-particle density matrix, Ak^k (t)- The qualitative behavior of Ak.k {t) is the same 
for the homogenous and the inhomogeneous coupling. Note that the different chain fillings for the same 
parameter set lead to a different overall scale for Ak^k- The comparison shows that the frequency and the 
size of the beating oscillations is smaller in the Hermitian time evolution than in the ratchet case. For an 
inhomogeneous hybridization, k is no longer a good quantum number. Apparently, the mode coupling in 
the Hermitian case noticeably reduces the fluctuations of Ak^k{t) around its time average. In the ratchet 
case, the Hamiltonian can still be written as a sum over commuting fc-modes, see eq. ( 14.241 ). Therefore, 
we observe a pronounced beating in Ak,k{t) even for inhomogeneous hybridizations. As we shall see in 
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Fig. 7 Diagonal part of the single-particle density matrix Ai^,k{t) for an inhomogeneous chain hopping, 7; = 
J exp(— Z /L), with L = 200, for the Hermitian time evolution (left) and the ratchet case (right). 

Sect. |5] interactions introduce a strong scattering between particles so that oscillations of the k modes are 
absent and Ak,k{t) quickly relaxes to its steady-state distribution. 

The non-diagonal components of Ak,p enter the expression for the current density ( I3.6l l. A small tran- 
sient current sets in at i > and decays to zero for long times Jt ^ 1. For our case study we choose 
7; — J exp{—l / L) so that we insert particles more readily in the left part of the chain than in the right. 
Originally, there is a cuiTent from the left to the right of the chain {jc > 0) which is reflected at the right 
chain end so that jc{t) changes sign as a function of time. As for the particle density, the current density 
shows oscillations which decay algebraically. 

In Fig. |8]we show the transient current density for the Hermitian time evolution and the ratchet case 
for L = 200 sites. The amplitudes differ for the same parameter set which is mostly due to the fact 
that there are more particles in the chain for the ratchet case. The frequency of the current oscillations 
is comparable. The current oscillations correspond to a swapping of the particles between the right and 
left chain boundaries. Since the mobile particles are concentrated around e(fco) — (fco = 7r/2), their 
dominant velocity is given by — {d£{k)) /{dk)\j.^^. — 2. Therefore, their travel time between the 
chain boundaries is T/2 w L/vq — L/2 so that the oscillation period is T w L. As seen in Fig. [8] 
this rule applies better for the ratchet case where the distribution Ak p{t) is more peaked around {k,p) ~ 
(7r/2,V2). 




Fig. 8 Current density jc{t) for non-interacting electrons for 7; = Jexp(— Z/L), unitary and non-unitary time 
evolution, for L = 200 sites. 
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4.2.2 Ratchet case 

Here we derive the formulae on which the data in Figs.|6]|2]and[8]are based. Even in the case of a non- 
symmetric coupling we may write the ratchet Hamiltonian in the form 



with 



Hk,Hp 



H = J2Hk , Hk^ EkCilak + al]3k (4.24) 

k=l 

for all fc, p. In eq. ( 14.241 ) we introduced the pseudo-fermionic operators 



Bk = 



7; sm 



1=1 



Trkl 
L + 1 



si , 



(4.25) 



which obey |i3fe, Bpj = js];, S+j = for all k,p. In particular, {Bkf = [Blf = 0. However, the 
operators and their adjoint do not anti-commute. 



7; I sm 



(=1 



nkl 
L + 1 



T:pl 
L + 1 



H{k,p) , 



(4.26) 



where we introduced the commutator function H{k,p). For ji — 7, we have Bk = 76^ and H{k,p) 
I7p4,p. 



The wave function can be obtained as in Sect. 14.21 

L 



\m) = l[{^ + 9k{t)Bkal)\m) , 9k{t) 



1 ~ e 



-iekt 



k=l 



£k 



(4.27) 



because the components Hk of H commute with each other. As shown in appendix|A] the single-particle 
density matrix can be written as 



■^p,q ~ 9q9p 



H{p,q)+ J2 Hip,ki)j:{k,,k2)H{k2,q) 
ki,k2—l 



(4.28) 



where E(/ci,A;2) are the entries of the self-energy matrix S. This is obtained from the solution of the 
matrix equation 



S = -(l + f •H)-^f . 



(4.29) 



Here, the entries of the matrices 1, H and f are la.b = 6a.b (unit matrix), Ha.b = H{a, b) (commutator 
matrix) and fa.b = Sa^bfa, fa = |.9aP (diagonal vertex matrix). The numerical solution requires the 
inversion of an X x L matrix. 



5 Interacting fermions 

In the thermodynamic limit and for non-interacting electrons, the average particle number in the chain goes 
to a constant at long-times only due to decoherence. On the other hand, a true relaxation of the particle 
density and other physical quantities can only take place when interactions are included. 

To investigate the model with the density-density interaction terms, eq. (13.31 1. included, we use an adap- 
tive time-dependent density matrix renormalization group (<-DMRG) algorithm. In DMRG the wave func- 
tion of the problem under consideration is approximated by a matrix-product state in a truncated Hilbert 
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space. The time evolution of this state is then calculated by using a second-order Trotter- Suzuki decom- 
position of the time-evolution operator. For details of the algorithm the reader is referred to Refs. f TSlflQl . 
In our calculations we use a Trotter slicing St = 0.1 and keep x — 250 states. It is known that the matrix 
dimension necessary to represent faithfully the time-evolved state increases exponentially with time. The 
i-DMRG algorithm is therefore only capable of describing the time evolution at short and intermediate 
times. 

5.1 Particle number in the chain 

In order to test the algorithm and the limitations for the time evolution, we perform numerical calculations 
for the particle density of non-interacting fermions with a homogeneous coupling where analytical results 
are available, see Sect. 14. II In Fig.|9]we compare numerical results obtained with the i-DMRG algorithm 
with exact results in the non-interacting case, both for the unitary and ratchet time evolution. As seen 
from the figure, the t-DMRG results (x — 250 states kept) are reliable for times Jt ^ 14 {Jt ^ 8) in the 
Hermitian (ratchet) time evolution. 




Fig. 9 Particle density in the chain nc{t) for non-interacting fermions on a chain with L — 50 sites (a) for the 
Hermitian time evolution (7/ J — 1) and (b) for the ratchet time evolution {-fin/ J = 1). Full lines; analytical result; 
dots: data from f-DMRG (x = 250 states kept). 

Apparently, the t-DMRG is not able to reach the limit of large times when the two-level systems will 
have dephased completely. However, the algebraic decay of the phase coherence in nc{t) as given in 
eq. ( 14.131 1 and as shown in Fig.|2]is special to non-interacting fermions. Interactions between the particles 
lead to a much faster dephasing and, moreover, to a true relaxation. In Fig.[TO]we show examples for the 
particle density in the chain for the Hermitian time evolution (7/./ — 1) in the presence of various types 
and strengths of the density-density interaction. In all cases, the interactions damp the oscillations and the 
average particle number reaches its limiting value much faster 

Depending on the type of interactions present, the energy of the system in its initial state 

Sini - (^-Oli^l^o) (5.1) 

is different. Since the Hamiltonian is a conserved quantity, the energy cannot change during time evolution. 
In the case where = (Fig. [TOh ) the system starts in a state with Ei^i = 0. Another state with zero 
potential energy for this type of density-density interactions is the twofold degenerate charge ordered state 
where chain and source sites are occupied alternatingly (Wigner lattice). This is shown pictorially in 
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Fig. 10 Particle density in the chain nc{t) for i = 50 lattice sites and 7/J = 1 for the Hermitian time evolution. 
We show t-DMRG results for non-interacting fermions (Vc = Vs = Vcs = 0, solid lines) in comparison with 
interacting fermions: (a) Vc = Vcs — 2,Vs — (dashed line) and Vc = Vcs = 8, K =0 (dot-dashed line); (b) 
Vc = Vcs = 0, K = 2 (dashed line) and Vc = Vcs =0,Vs=8 (dot-dashed line). 



Fig. [TT] At least in the limit where Vc = Vcs is large, i.e., where the potential energy dominates over 
the kinetic energy, we therefore expect the system to reach a stationary state which is close to this charge 
ordered state. Therefore, the average particle density approaches nc{Jt ^ 1) = 1/2 (see Fig. [TOk) for 
large interactions which is larger than in the non-interacting case. Apparently, the Coulomb interaction 
reduces the degree of delocalization of electrons on the chain whereby the Pauli blocking becomes less 
effective. The reduction of the Pauli pressure increases the density of fermions in the chain. 

The case Vc — Vcs — with Vs (Fig. [TOb) has an energy Eini — {L — l)Vs. For the situation 
where Vs is smaller than the bandwidth, energy conservation can be fulfilled by compensating the loss of 
potential energy by kinetic energy which the fermions gain by entering the chain. Numerically we find 
that the density of these almost free fermions in the chain can reach ndJt ^ 1) sa 1/2, see Fig.fTOb. For 
Vc = Vcs ~ with Vs larger than the bandwidth, on the other hand, there is no possibiUty of compensating 
the loss in potential energy with kinetic energy in the chain. Hence the particles are almost all confined to 
the bath sites by energy conservation and nc{Jt ^ 1) remains small, see Fig.fTOb. 



o 



o 



o 



-o- 



-o- 



-o- 



Fig. 11 Charge ordering in a Wigner crystal state. The crosses mark the occupied sites. 



The fast damping of the oscillations of the particle density is a genuine interaction effect which does 
not depend on the formation of a charge-density wave. This can be seen in Fig. [12] where we show the 
time evolution of the particle density in the chain for Vs = Vc = Vcs = 1,2,3,4 and 7 = J = 1. As 
for Vc — Vcs — 0, Vs ^ 0, due to the large initial energy, E^ini ^ {L — l)Vs, of the system the number 
of accessible states is severely restricted. Therefore larger interactions, V > 3, not only lead to a faster 
relaxation but also suppress the particle density in the chain. 
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Jt 

Fig. 12 Particle density in the chain nc{t) for L — 50 lattice sites and 7/ J = 1 for the Hermitian time evolution. We 
show t-DMRG results for interacting fermions: Vc ~ Vs = Vcs = 1 (solid line), Vc ~ Vcs = Vs — 2 (dotted line), 
Vc = Vcs = Vs = 3 (dashed line), and Vc = Vcs = 14 = 4 (dot-dashed line). 

In the limit where Vg — Vc ^ Vcs is much larger than the bandwidth, only the two states with the 
largest potential energies are reachable due to energy conservation: a full bath or a full chain. In general, 
the system can tunnel between these two possible states. For a large system size, however, the probability 
of all particles coherently tunneling is very small and therefore the associated timescale is very long. If 
we reduce the system size to L = 2, on the other hand, the oscillations between these states occur on a 
relatively short timescale as can be seen in Fig.[T3] For L — 2 the Hamiltonian can be exactly diagonalized 
for arbitrary interactions and we show here exact results. 




Fig. 13 Particle density in the chain nc{t) for L = 2 lattice sites and 7/J = 1 for the Hermitian time evolution. We 
show exact results for non-interacting and interacting fermions: Vc = Vs = Vcs ~ (solid line), Vc = Vcs = 14 = 8 
(dashed line). 

Next, we consider the ratchet case with interactions. Note that time evolution is limited to Jt < 6 so 
that we cannot follow the time evolution for too long. In Fig. [14] we show t-DMRG results for the ratchet 
case with jin/ J — 1- Again, the interactions damp the fluctuations and, for the specific cases considered 
here, lead to an increase of the average particle number in the chain. 

5.2 Single -particle density matrix 

We will concentrate on the Hermitian time evolution for which i-DMRG data are available for short and 
intermediate times. In the interacting case, scattering processes between fermions provide an efficient 
mechanism for the relaxation of the single-particle density matrix. Therefore, we expect that Ak.k{Jt ^ 1) 
becomes independent of time. Note, however, that this distribution may, but need not, correspond to a 
thermal distribution for an appropriately chosen ensemble at a temperature determined by (if) thermal = 
(if) initial. The interesting question whether or not this thermal state is reached will be addressed in a 
forthcoming publication Il20l . 

In Fig. [15] we show results for interaction parameters Vc — Vcs = 2, T4 = and Vc — Vcs — 0, 
Vg = 2. For Vc = Vcs = 2, T4j = at short times, the single-particle density matrix is still very close 
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Fig. 14 Particle density in the chain nc{t) for L = 50 lattice sites and fin/ J = 1 for the ratchet time evolution. 
We show t-DMRG results for non-interacting fermions (Vc ~ Vs — Vcs ~ 0, solid lines) in comparison with 
interacting fermions; (a) Vc = Vcs = 2,Vs =0 (dashed line) and Vc = Vcs = 4, Vs =0 (dot-dashed line); (b) 
Vc = Vcs = 0, Vs = 2 (dashed line) and Vc = Vcs = 0, Vs = 4 (dot-dashed line). 
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Fig. 15 Single-particle density matrix Ak,kit) for the Hermitian time evolution with 7/,/ — 1 for a chain of length 
L = 50 with interaction parameters Vc = Vcs = 2, Vs = (left) and Vc — Vcs = 0, Vs = 2 (right). 



to the non-interacting case shown in Fig. H) At longer times, the interactions lead to a rapid damping of 
the oscillations and, already at J< = 10, the distribution in both cases becomes almost time-independent. 
For Vc = Vcs = 2, Vs = the distribution is very flat Ak^k{t ^ 1) — 1/2- This observation supports 
our interpretation of the data shown in Fig. [TOh : for strong interactions, the fermions tend to form a 
Wigner lattice of alternatingly occupied chain and source sites which corresponds to a completely flat 
single-particle density matrix. For Vc — Vcs — 0, Vs = 2, on the other hand, the distribution at long 
times resembles more a thermal Fermi distribution. Here the loss of potential energy due to particles 
leaving the bath has to be compensated by a gain in kinetic energy. For the interaction strength considered, 
approximately half of the particles enter the chain in the long-time limit, see Fig.fTOb. These particles are 
almost free thus explaining why the long-time distribution shown in Fig. [15] is qualitatively similar to a 
thermalized free fermion distribution. 
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In contrast to the non-interacting case shown in Fig. |4] Ak,k{t) is a smooth function of the pseudo- 
momentum k. Despite the fact that our t-DMRG cannot go much beyond Jt — 10, our data suggest that 
the electron-electron interactions lead to a fast relaxation of the single-particle density matrix to its value 
in the long-time limit. 

5.3 Current density 

For non-interacting particles, the current decays to zero algebraically in time. The reason is the dephasing 
of the fc-modes. In the presence of interactions, the current not only dephases but decays (exponentially) 
quickly. In Fig.[T6]we show the result for the Hermitian time evolution for L = 50 sites and two sets of 
interaction parameters. The fast decay of the current is clearly observed. 




Fig. 16 Current density jc(i) for interacting electrons for 7; — Jexp(—Z/L) in the Hermitian time evolution (J = 1) 
for L = 50 sites with Vc = Vcs = 2, Vs = (circles) and Vc ~ Vcs — 4,Va = (squares). For comparison we also 
show the result for noninteracting fermions (full line). 

As discussed in Sect. 14.21 the inhomogeneous injection of non-interacting fermions results in a swapping 
of particles between the chain boundaries. As seen in Fig. [8] the cuiTent density increases to a maximal 
value of Jcinax'"* ~ 0.1 at w 50, oscillates with period T w i, and decays to zero algebraically 
slowly. For interacting fermions, however, the current density reaches its maximum of Jc"niax ~ 02 at 
Jt « 5, and rapidly decays to very small values already at Jt « 10. Apparently, the interactions prevent 
the buildup of a coherent mesoscopic particle wave package, and only local transient current fluctuations 
can be seen. Note that this effect is most prominent in one dimension where particles can scatter forward 
and backward only. In addition, we have chosen quite sizable interactions to make the rapid relaxation 
visible within our restricted time interval in t-DMRG. For a more detailed analysis, longer times must be 
studied. This is beyond our numerical capabilities. 

6 Summary and Conclusions 

In this work, we studied the injection of spinless fermions from filled reservoir sites into an empty chain 
on which the particles can move between neighboring sites. We treated the cases of homogeneous and 
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inhomogeneous hybridizations between the source and chain sites under a Hermitian time evolution and 
for the quantum ratchet case where the chain particles cannot return to the source sites. For non-interacting 
electrons, we provide analytical expressions as a function of time for the particle density, the single-particle 
density matrix, and the current density. In the presence of interactions, we use the t-DMRG algorithm to 
describe the time evolution of observables for short to medium time scales. 

The Hermitian time evolution and the ratchet time evolution both show some qualitatively similar effects 
which is quite surprising at first. Despite the fact that the particles cannot return from the chain to the source 
sites in the ratchet case, their injection remains a coherent quantum-mechanical process. Therefore, the 
well-known quantum-mechanical phenomena such as constructive and destructive interference and Pauli 
blocking are seen in our observables. In particular, there are in both cases oscillations around the long- 
time average of the chain filling. Naturally, the results differ quantitatively. For example, the chain filling 
(transient current) is larger (smaller) for the ratchet case than for the Hermitian case for the same parameter 
set. 

The time evolution of non-interacting particles is characterized by the absence of relaxation so that 
fairly large dephasing times, Jt > 10'^, dominate the time evolution of our quantities. For large systems, 
the dephasing ultimately leads to steady-state expressions for the particle and current densities but the 
single-particle density matrix reveals that the spectrum is that of a single-particle Hamiltonian. For a ho- 
mogeneous coupling, the system consists of independent two-level systems and the single-particle density 
matrix reflects the beating of these two-level systems. For inhomogeneous hybridizations, the beating of 
the diagonal terms Ak,k{t) is less pronounced because k is no longer a good quantum number The inho- 
mogeneous coupling between chain and source sites generates non-diagonal elements in the single-particle 
density matrix which, however, are small and dephase in time. They give rise to the coherent motion of 
a spreading particle wave packet between the chain ends which results in a transient oscillating current 
density. 

When we introduce substantial interactions between the particles, the picture changes qualitatively. 
Depending on the type of density-density interactions considered, the system starts with different initial 
energies. Since the energy is conserved during time evolution, only parts of the Hilbert space are acces- 
sible. Interactions can therefore increase or decrease the particle density in the chain compared to the 
non-interacting case. Furthermore, the single-particle density matrix in the long-time limit can show a 
completely flat distribution in cases where the potential energy dominates, or become similar to a thermal- 
ized free fermion distribution when the kinetic energy plays the dominant role. Most importantly, instead of 
a slow dephasing or beating phenomena, we observe a rapid relaxation to stationary values on time scales 
Jt = 10, not only for the particle and current densities but also for the single-particle density matrix. This 
shows the importance of interactions for a fast relaxation for all physical quantities of interest. In how far 
the obtained time-independent results do coincide with those predicted by thermodynamic considerations 
is a question we are planning to address in more detail elsewhere [20 1. 

Acknowledgements J.S. acknowledges support by the MAINZ (MATCOR) school of excellence and the DFG via 
the SFB/Transregio 49. 

A Single-particle density matrix for non-interacting fermions (inhomoge- 
neous ratchet case) 

A.l State 

As seen from eq. ( 14.27b . the time-dependent state is given by 



L 




(A.l) 
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where \ip{0)) = Y\ sl\vac) and gk{t) ^ {1 — e '^''*)/efc. We rewrite the product as 

1=1 



\m) 



k 

ki<k2 

■ 51 9k,9k29k3Bk,al^Bk2al^Bk3ai^ 

ki <k2<k3 

\\m)- 



(A.2a) 
(A.2b) 

(A.2c) 

(A.2d) 



The first term, eq. (IA.2al i. describes a state where there is no particle on the chain, the second term, 
eq. ( IA.2bb describes the situation with one particle on the chain, and so on. 



A. 2 Numerator 

In the density matrix g, the operator aJ^Og does not change the number of particles on the chain. There- 
fore, 



m)\ald^\'4>{t)) = 9q9;{BlBg)o 

+ 9lM^k2BlakMMlBA)^ 



(A.3a) 
(A.3b) 



9q9k,9k2 



E 



ki<k2^{q,p) 



D Si,S2,S3 



S3 



(A.3c) 

Sl^S2i^S3^Si 

X (a^3 a,, Bl^ as, B]^ Bk, i?fe2 ^i^l^ 



We can factor out the expectation value for the chain operators and use Wick's theorem flT] so that we can 
identify some of the momenta. Relabeling some indices leads to the following expressions for the second 
and third term of the numerator. 



^AM = J2 9q9k9*k9;{BlBlBkB,)o , 
k=i'={p,q) 

^E3^ = 9q9;\9kA'\9k2\'{BlBlBlBk,Bk2B,)o 



E 



(A.4) 
(A.5) 



The expectation values are calculated using Wick's theorem. To this end we need {i3jji3q)o — H{p, q), see 
eq. ( 14.26b . The expectation values in eqs. ( IA.4b and (IA.5I ) can be expressed as determinants. 



{BlBlBkB. 



q)0 



H{p,q) Hip,k) 
H{k,q) H{k,k) 



(A.6) 



{BlBl2BlBk^Bk2B,), = 



H{p,q) H{pM) H{p,k2) 
H{ki,q) H{ki,ki) Hik^M) ■ (A.7) 
H{k2,q) H{k2,ki) H{k2M) 
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This leads to the result 

m)\aia^\m) 



9q9\ 



H{p,q) H{p,k) 
H{k,q) H{k,k) 



(A.8) 



fei ,k2 



Hip,q) H{pM) HipM) 
H{ki,q) H{kiM) Hiki,k2) 

H{k2,q) H{k2,ki) ff(fc2,fc2) 



We were allowed to drop the restriction on the summations due to the properties of the determinants. 
A. 3 Self-energy 

As in all perturbative calculations, we can translate the terms in (IA.8b into diagrams 1211 with 'inter- 
nal vertices' k of strength fk = \gk\'^ and lines between vertices ki and fc2 representing the factors 
H{ki, fc2). In the thermodynamic limit, L — > oo, the linked-cluster theorem applies (2T\, i.e., the de- 
nominator {ip{t)\ip{t)) in Ap^q cancels the unconnected diagrams in the numerator In Ap^q we have to sum 
over aU diagrams which connect the external vertex p with the external vertex q. We find 



A * 

P. 9 ~ 9q9p 



i7(p,(Z)-^i7(p,fc)/feF(fc,g)+ Hip,k,)fk,H{ki,k2)fk2H{k2,q)T ■ 



ki ,k2 



We introduce the self-energy Xl(a, b) il2T1l and prove eq. < l4.28ll . 



(A.9) 



A 



p,q — 9q9p 



Hip,q)+ J2 H{pM)nkiM)H{k2,q) 



(A. 10) 



ki,k2 



where the self-energy entry S](a, h) is obtained by summing all possible paths from the internal vertex a to 
the internal vertex b with altogether n internal vertices {n— 1,2,...) and (n — 1) lines between them. 



E(a,6) 



^a,bfa ~t~ fa 



H{a, b)h - H{a, k)hH(k, b)h 



Y H{a,ki)fk,H{ki,k2)fk2H{k2,b)fbT 



ki,k2 



(A. 11) 



The solution of this matrix equation provides S, see eq. ( 14.291 ). 
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